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Abstract. Structural and energy characteristics of the smallest magnetic endofullerene Fe@C2o 
were calculated using the density functional theory. The ground state of Fe@C2o was found to 
be a septet state, and the magnetic moment of Fe@C2o was estimated to be 8 Bohr magnetons. 
The characteristics of an (8,8) carbon nanotube with a single Fe@C2o inside were studied with a 
semiempirical approach. The scheme of a magnetic nanorelay based on cantilevered nanotubes 
filled with magnetic endofullerenes was examined. This nanorelay is turned on as a result of 
bending of nanotubes by a magnetic force. The operational characteristics of such a nanorelay 
based on (8,8) and (21,21) nanotubes fully filled with Fe@C2o were estimated and compared to 
the ones of a nanorelay made of a (21,21) nanotube fully filled with experimentally observed 
(Ho3N)@Cgo with the magnetic moment of 21 Bohr magnetons. The room-temperature opera- 
tion of (21,21) nanotube-based nanorelay s was demonstrated. 

Keywords: single-walled carbon nanotubes, endofullerenes, magnetically operated nanorelay 
1 INTRODUCTION 

The remarkable elastic properties and metallic conductivity of carbon nanotubes (CNTs) allow 
applications of nanotubes as parts of nanoelectromechanical systems (see, e.g., Refs. 1-3). In 
fact, a wide set of electromechanical nanorelays based on the crossbar relative position of two 
carbon nanotubes [4], cantilevered [5-8] and suspended [9] nanotubes, telescopic extension of 
nanotubes [10] and a magnetic shuttle inside a carbon nanotube [11] has been implemented. 
Operational principles of these nanorelays and possibility to use them as memory cells are 
considered in Ref. 12. 

A variety of endofullerenes and nanotubes filled with fullerenes (nanotube peapods) are 
obtained in macroscopic amounts [13], including nanotube peapods filled with magnetic endo- 
fullerenes [14-17]. In the present work we propose a new type of a nanorelay based on nanotube 
peapods encapsulating magnetic endofullerenes. Operation of this new nanorelay is determined 
by the balance of three major forces: magnetic, elastostatic, and van der Waals ones. 
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Endofullerenes with encapsulated iron atoms were observed experimentally [18, 19]. Re- 
cently the magnetic moment of the endofullerene Fe3@C6o was calculated [20]. Here we con- 
sider the fullerene C20 with encapsulated iron atom as the simplest model of a magnetic endo- 
fullerene. We show the possibility that the smallest magnetic endofullerene Fe@C2o can exist. 
The density functional theory (DFT) approach is used here to study the structure, energetics 
and the magnetic moment of the endofullerene Fe@C2o- The possibility of encapsulation of 
small fullerenes inside nanotubes was studied in Ref. 21. It was shown that an armchair (8,8) 
nanotube is the smallest carbon nanotube (CNT) which can encapsulate fullerene C20 [21]. In 
the present work the binding energy, the equilibrium position, and energetic barriers for mo- 
tion and rotation of the endofullerene Fe@C2o inside the (8,8) nanotube are calculated using a 
semiempirical potential. 

The magnetic moment was measured for a wide set of fullerenes with encapsulated magnetic 
atoms [17,22-27]. It was also shown that the magnetic moment is larger for magnetic endo- 
fullerenes inside carbon nanotubes than for the same isolated magnetic endofullerenes [17]. 
Nowadays the largest value of the magnetic moment of 21 Bohr magnetons was observed for 
(Ho3N)@C8o [24]. Here operational characteristics of the proposed magnetic nanorelay are 
calculated for the cases of nanotubes filled with the smallest theoretically possible magnetic 
endofullerene Fe@C2o and the magnetic endofullerene (Ho3N)@C,so with the largest observed 
magnetic moment. 

The paper is organized as follows. Section 2 presents the DFT calculations of properties of 
the endofullerene Fe@C2o- Section 3 is devoted to characteristics of endofullerene encapsula- 
tion inside the (8,8) carbon nanotube. Section 4 contains the principal scheme of the nanotube 
peapod-based nanorelay controlled by magnetic field and estimations of its operational charac- 
teristics. Our conclusions are summarized in Sec. 5. 

2 ENDOFULLERENE Fe@C 20 

The study of carbon-iron systems is important for physics and chemistry of carbon nanostruc- 
tures. It is known that iron is a catalyst for synthesis of carbon nanostructures and diamond. 
Thus, in the process of CNT synthesis [28, 29] iron nanoparticles were sometimes left on free 
tips of nanotubes. The mechanisms of such processes are still not established. We consider 
behavior of an iron atom in a quantum-size structure, fullerene C20, as an example of one of the 
simplest carbon-iron systems. 

According to the quantum chemical calculations [30], C atoms of the fullerene C20 are lo- 
cated at 0.205 nm from its center (for the average C-C bond length of 0. 146 nm). The iron atom 
in the ferrocene CioHioFe is located at the same distance from the carbon atoms (see, e.g., [31]). 
Due to this fact we assume that an iron atom can be put inside the fullerene C20. The iron atom 
in the center of the fullerene C20 is located in the force field with the approximate symmetry 
of icosahedron J/j. The field with such symmetry does not split the D terms of atoms [32, 33], 
therefore the iron atom can preserve its angular momentum inside the fullerene C2o- 

Spin states of the endofullerene Fe@C2o have been calculated using the spin-polarized den- 
sity functional theory implemented in the NWChem 4.5 code [34] with the Becke-Lee- Yang- 
Parr exchange-correlation functional (B3LYP) [35, 36]. 18 inner electrons of the iron atom are 
emulated with the help of the effective core potential — CRENBS ECP [37] (only 8 valence s-d 
electrons were taken into account explicitly). The 6-3 1G* basis set is used to describe electrons 
of the carbon atoms. 

The calculated energies of the quintet (5 = 2), triplet (5* = 1) and singlet (5 = 0) states of 
the endofullerene Fe@C2o are found to be 0.61, 0.86 and 1.21 eV higher than the energy of the 
ground state. The ground state of the endofullerene Fe@C2o is found to be a septet state (with 
the total spin moment 5 = 3 and the multiplicity 25+1 = 7) and has the C211 symmetry. Since 
the ground state of a free iron atom [38] has the total moment J = 4 and the multiplicity 
25+1 = 5, the atomic magnetic moment of the iron atom inside the fullerene C20 brings about 
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Fig. 1. The calculated structure and spin density of atoms in units of h/2 of the ground state of 
the endofullerene Fe@C2o. The lengths of bonds (in angstroms) are in italics. The dotted line 
C*2 represents the two-fold symmetry axis. 

the ferromagnetic ordering of 7r-electrons of the unfilled electron shell of the fullerene. The 
magnetic moment of an iron atom is M Fe = —gj/j, s , where /i B = 5.788- meV/mT is the 
Bohr magneton and g = 1 + [ J( J + 1 ) - L(L + 1 ) + S (S + 1 )] / [2 J( J + 1 )] is the Lande splitting 
factor (spectroscopic splitting factor), J is the total angular momentum of an iron atom. For the 
ground state 5 £>4 term of the iron atom we use [38]: J = 4, L = 2, 5 = 2 and the Lande 
factor g = 3/2, this yields Mp e = gyj J {J + 1)^b ~ 6/ib- We assume that the endofullerene 
Fe@C2o magnetic moment is made up of the iron atom magnetic moment and the spin magnetic 
moment of the carbon cage (M cc ps 2/xb). So we obtain M e { w Mve + A/ cc = 8/jb- 

The calculated structure of the ground state of the endofullerene Fe@C2o is shown in Fig. 1. 
The iron atom is located in the center of the endofullerene. The distances between the iron atom 
and the carbon atoms are 0.206 and 0.210 nm for the carbon atoms shown in Fig. 1 by the open 
and filled circles, respectively. Note that these distances are smaller for the carbon atoms with 
anti-parallel to the iron magnetic moment spin density (attracting to the iron atom) and larger 
for the carbon atoms with parallel spin density (repulsing from the iron atom). 

The Mulliken population analysis [39] is used to calculate the charge distribution of Fe@C2o- 
It is found that valence electrons of the iron atom partially shift to the carbon cage, and the iron 
atom gains the charge of 1 .2e, while the carbon atoms located closer to and farther from the iron 
atom gain the charges of — 0.04e and — 0.07e, respectively, where e is the elementary charge. 

Our calculations of the empty fullerene C20 show that it is a fairly stable system with the 
binding energy of 6.15 eV/atom (slightly less than 7.0 eV/atom [40] for the fullerene C60) and 
the singlet ground state (the triplet state is 65 meV higher). Though, isomers of cluster C20 
with lower energy may exist, calculations [41] show the high thermal stability of the fullerene 
C2o- Up to now the fullerene C20 was detected only in the gas phase [42]. According to our 
calculations the incorporation of an iron atom into the fullerene C20 requires 5.1 eV. This is 
less than the C atom binding energy in the empty fullerene C20, which indicates the possibility 
of Fe@C2o existence, while the stability of such a system remains at the same level as for the 
empty fullerene C2o- 
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3 Fe@C 20 INSIDE CNT (8,8) 

The structure of the endofullerene Fe@C2o described in Sec. 2 is used to find the ground state 
and to study motion of this endofullerene inside the (8,8) carbon nanotube with all equal bonds 
0. 1423 nm in length. The van der Waals interaction between carbon atoms of the endofullerene 
and the nanotube at the interatomic distance r is described by the Lennard- Jones 12-6 potential 
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with the parameters e = 2.755 meV and a = 0.3452 nm. These parameters of the Lennard- 
Jones potential for the fullerene-nanotube interaction are obtained as the average values of the 
parameters [43] for fullerene-fullerene and fullerene-graphene interactions, in accordance with 
the procedure described in [43]. The cut-off distance r = r c of the Lennard- Jones potential 
is taken equal to r c = 1.5 nm. For this cut-off distance the errors of our calculation in the 
interaction energy ?7w between the endofullerene Fe@C2o and the (8,8) nanotube and in the 
barriers for relative motion and rotation of the endofullerene inside the nanotube are less than 
0.1%. Both the fullerene and the nanotube are considered to be rigid. Influence of structure de- 
formation is not essential both for the interwall interaction of carbon nanotubes [44,45] and the 
intershell interaction of carbon nanoparticles [46,47]. For example, taking into consideration 
the structure deformation of the shells of the C6o@C24o nanoparticle only slight changes of the 
barriers for relative rotation of the shells are found (less than 1%) [46,47]. 

The ground state interaction energy between the endofullerene Fe@C2o and the (8,8) nano- 
tube is calculated to be —1.598 eV. The position of the endofullerene Fe@C2o inside the (8,8) 
nanotube at the ground state is shown in Fig. 2, displaying the angle of 52° between the sym- 
metry axes of the endofullerene and the nanotube. The contour map of the interaction energy 
U-w(<p, z) is presented in Fig. 3, where (p is the angle of relative rotation of Fe@C2o about the 
nanotube axis, and z is the relative displacement of Fe@C2o along this axis. Figure 3 shows 
that both barriers for motion of the endofullerene along the nanotube axis and for rotation of 
endofullerene about this axis are 0.64 meV. The minimal barrier for motion of endofullerene 
along the nanotube axis corresponds to the case where this motion is accompanied by rotation 




Nanotube axis 



Fig. 2. The calculated structure of the ground state of the (Fe@C2o)@CNT(8,8) peapod nano- 
tube. The diameters of the endofullerene and the nanotube are indicated. Thinner lines corre- 
spond to the back side. 
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Fig. 3. The contour map of the interaction energy Uy/(<p, z) between the endofullerene Fe@C2o 
and the (8,8) carbon nanotube as a function of the relative displacement z (in nm) of Fe@C2o 
along the nanotube axis and the angle ip (in degrees) of relative rotation of Fe@C2o about it. The 
equipotential lines are drawn at an interval of 0.05 meV. The maximum f/w value is 0.64 meV. 

about the nanotube axis. Figure 3 shows that the value of this minimal barrier between adjacent 
minima is 0.45 meV. 

The calculations performed here show that free motion of the endofullerene along the nano- 
tube axis occurs at room temperature. Therefore a single endofullerene Fe@C2o is able to move 
to a nanotube end where it can be trapped in the potential well associated with the interaction 
between the endofullerene and the nanotube cap. In much the same way endofullerenes are able 
to assemble in a one-dimensional cluster inside the nanotube due to van der Waals interaction 
between the endofullerenes. 



4 NANOTUBE-BASED MAGNETICALLY OPERATED NANORELAY 

We propose a new type of a nanorelay based on peapod nanotubes with encapsulated magnetic 
endofullerenes (see Fig. 4) which allows to disconnect control and controlled circuits, analo- 
gously to known [48] macroscopic reed switches. The proposed nanorelay has the following 
operational principles. When a magnetic field is applied, the majority of the magnetic moments 
of the endofullerenes lines up in the direction of the magnetic field as it is shown in Fig. 4(a). 
As a result of the magnetic moment alignment, the attraction between the nanotubes arises. The 
condition for the nanorelay switching is determined by a balance of the magnetic and elasto- 
static forces. 

The interaction energy of two magnetic dipoles with the moments Mj and M 3 - of the endo- 
fullerenes encapsulated in the right and left nanotubes, respectively, is given by (see, e.g., 
Ref. 49) 

(MjMj)-3(nM^ 
Uij — Mo ^3 j (A> 

where yu = 47r-10~ 7 H/m is the magnetic vacuum permeability, r is the distance between the 
endofullerenes, and n is the unit vector parallel to the line joining the centers of the endofulle- 
renes. 
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Fig. 4. The scheme of the magnetic nanorelay based on nanotubes: 1 is the controlled circuit, 
2 is the electrode with the attached peapod nanotube, 3 is the wire passing the current I through 
it and used for inducing the magnetic field with the induction B, which brings tips of nanotubes 
into contact with each other. 

4.1 Single magnetic endofullerene at the tip of nanotube 

Let us consider the case where each nanotube of the nanorelay contains a single magnetic endo- 
fullerene Fe@C2o at the tip of the (8,8) nanotube. In this case the magnetic force F m of inter- 
action between the nanotubes is given by the result of differentiation of Eq. (2) with respect to 
2x. This gives 
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where f m is the force of interaction between two magnetic moments Mj 
which line up in the direction of the magnetic field and N m < 1 is the fraction of the total 
number N e f = 1 of the endofullerenes with the magnetic moments parallel to the magnetic 
field. This number N m is determined by the Langevin function C(M e f) and depends on the 
magnetic field induction B and temperature T (see, e.g., Refs. 49 and 50): 



N m = N*L{M*) = N ef 



coth 



( MgB 

\ fc B T 



M ef B 
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where N e f = 1 is the total number of the endofullerenes inside each nanotube and k s is the 
Boltzmann constant. Since the considered magnetic moments are large we use the Langevin 
function C instead of the Brillouin one. 

The elastostatic force acting on the carbon nanotube is given by 



k(x D — x), for the right nanotube, 
k(x — x ), for the left nanotube, 



(5) 



where x is the deflection of the nanotube tip, 2x a is the distance between the axes of the nano- 
tubes in the nanorelay in the absence of a load and k is the bending stiffness coefficient. The 
macroscopic approach is conventionally used to describe the elastic behavior of nanorelays 
based on cantilevered nanotubes (see, e.g., Refs. 51-53). For a load applied to the nanotube tip 
the bending stiffness coefficient is given by [54] 



2YI a 
L 3 



(6) 



where Y is Young's modulus, determined to be approximately 1.2 TPa [55-57], L is the nano- 
tube length, I nt = (7r/4)(i?£ X — Rf n ) is the moment of inertia of the nanotube cross-section, and 
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Fig. 5. The schematic representation of the dependencies of the magnetic force F m and the 
elastostatic force F s on the coordinate x of the nanotube tip, 2x a is the distance between the 
nanotube tips when magnetic field is absent. The nanorelay is turned on if F s < F m at any coor- 
dinate x, i.e., the distance between the nanotube tips should be less than 2x for the nanorelay 
turning on by the magnetic field. 

i?i n and R ex are the inner and outer radii, respectively, of the nanotube. The moment of inertia 
of a single-walled nanotube can be estimated as 

/„, = ttR 3 AR, (7) 

where R is the nanotube radius and Ai? = i? ex — i?; n sa 0.34 nm is the effective thickness of the 
single-walled nanotube, which is taken to be equal to the interlayer distance of graphite; i? ex = 
R + 0.17 nm and i?; n = R — 0.17 nm. The radius of an (n, m) nanotube is determined by the 
following expression [58,59]: R = (ac^c / 2-k) 2>(n 2 + to 2 + nm), where clq-c = 0.142 nm 
is the C-C bond length in the nanotube. For the (8,8) nanotube R = 0.542 nm. 

The nanorelay is turned on by the contact of nanotube tips. The nanorelay turning on is 
possible if the magnetic force is larger than the elastostatic force at any deflection x of the 
nanotube tip. Thus the equalities of magnetic and elastostatic forces F % (x{) = F m (x{) and 
their derivatives dF s /dx\ x=Xi = dF m /dx\ x=Xi determine the maximum value of the distance 
2x between the nanotubes (for defined iV m ) for which the nanorelay turning on is possible 
(see Fig. 5). Substituting Eqs. (3) and (5) into these equalities, we obtain the following set of 
equations 

2 xf 8 x\ 

4 Ar 32 xl /2 Vmk 

Xi = -x , N m = —— — , (9) 

5 375 M e f y/flo 

where Xj is the coordinate of the nanotube tip when F m = F s and x a <C L. 

The amplitude of thermal oscillations of the nanotubes has to be small enough to preclude 
the possibility of spontaneous closure of the nanorelay in the absence of the magnetic field. The 
nanorelay cannot operate in the case of the average thermal oscillation amplitude larger than 

— £™ n , where 2x™ n = 2R+ Ai? is the distance between the nanotubes which corresponds to 
their contact. For this condition one can use an estimate fc(x ,t h — x-™ n ) 2 /2 > k^T/2, where 2a;* 
is the minimal distance between the nanotubes limited by thermal oscillations. More accurately, 
by taking the expression for the average thermal oscillation amplitude [60], this condition is 
given by 

(x*-xrf> 0.4243^^, (10) 



Hence it follows that 



where k^T is thermal energy. 
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Fig. 6. (a) The dependence of the minimal magnetic field induction £> m j n on the nanotube 
length L for the nanorelay operated at T = 0. 1 K and based on the (8,8) carbon nanotubes (left 
axis). Each nanotube contains single magnetic endofullerene Fe@C2o with M e f = 8^b at the 
tip. The dependences of half of the minimal (2a;*) and maximal (2x™ ax ) distances between the 
nanotubes for which the nanorelay operation is possible on the nanotube length L (right axis); 
x mm _ 0.712 nm. (b) The dependence of the induction B necessary for the nanorelay operation 
with respect to the distance 2x a between the nanotubes for the nanotube length L = 0.5 fim. 
Arrows indicate for this L the values of x*, x™ ax and B^ shown in Fig. 6(a). 



We obtain from Eq. (9) at N m = 1 the maximal distance 2x™ ax between the nanotubes for 
the nanorelay turning on for a given nanotube lengths: 



V 32 

Substituting Eq. (9) and x = a;* (according to Eq. (10)) into Eq. (4), the minimal magnetic 
field induction B m \ n which is necessary for the nanorelay turning on can be numerically calcu- 
lated. The dependence of the minimal magnetic field induction B m [ n on the nanotube length L 
for the nanorelay based on the (8,8) nanotubes operated at T = 0. 1 K is shown in Fig. 6(a). Fig- 
ure 6(a) also shows the dependence of half of the minimal (2a;*) and maximal (2x™ ax ) distances 
(according to Eqs. (10) and (11), respectively) between the nanotubes for which the nanorelay 
operation is possible with respect to the nanotube length L corresponding to the case under con- 
sideration. Figure 6(b) shows the dependence of the induction B necessary for the nanorelay 
operation on the distance 2x Q between the nanotubes for nanotube length L = 0.5 /zm. 
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The measurements of van der Waals forces between the caps of the nanotubes about 1 nm 
in diameter gave the maximal value ps 1-2 nN [61]. For the considered case of the interaction 
between the lateral surfaces of nearly parallel nanotubes the van der Waals forces should be 
significantly larger. Whereas for the (8,8) nanotube 0.5 pm in length the bending stiffness 
coefficient is equal to k = 5 nN/mm, and the elastostatic force for the deflection of the nanotube 
tip, x = 1.2 nm, is equal to F s « kx = 6-10 -6 nN. Thus the van der Waals forces are 
significantly larger than the elastostatic forces. The nanorelay turning off can be achieved by 
applying an electrostatic potential (without current) to the same wire used for the nanorelay 
turning on [7]. 

The switching time r of the nanorelay can not be considerably less than the period of free 
bending vibrations of the nanotube. The frequencies of the free bending vibrations of the nano- 
tube are given by [60] : 

u '- 2 *f>-2V)l p ' (12) 

where j3 s for a given mode s is found from the solution to the equation cos /3 S cosh f} a + 1 =0; 
/So ~ 1-8751 for the fundamental mode, p m 2.15 g/cm 3 is the density of the material in 
the nanotube wall. For the (8,8) nanotube (R ex = 0.712 nm, i?; n = 0.372 nm) of the length 
L = 0.5 pm we estimate: 

2 



1 2L 

TW r 2 ^f(c + ^r 50ns ' (13) 

The induction B of the magnetic field induced by the current I at the distance d from the wire 
is [62]: 

B=^L. (14) 

To turn on the nanorelay by applying the magnetic field with the induction of 60 mT (for the 
(8,8) nanotube of 0.5 pm in length, see Fig. 6) it is necessary to use the current I = 300 mA 
passing through the wire 3 (see Fig. 4), which is positioned at distance d = 1 pm from nanotu- 
bes, for the time interval somewhat longer than r w 50 ns. 

4.2 Nanotube fully filled with magnetic endofullerenes 

Let us consider the magnetic interaction between two parallel peapod nanotubes of the length L 
fully filled with magnetic endofullerenes Fe@C2o- The number N m 3> 1 of the endofullerenes 
with magnetic moments parallel to the magnetic field is determined by the total number of 
encapsulated magnetic endofullerenes N e f 3> 1, the magnetic field induction and temperature. 
By summing up the forces = VE/ij, where J7y is given by Eq. (2), over all pairs of such 
magnetic dipoles, it can be easily shown that for the case where the distance 2x between the 
nanotubes is considerably less than the nanotube length L, 2x <C L (i.e., neglecting the tube 
edge effects), the magnetic force of interaction between the nanotubes takes the form 



Nm N m 



dUij ^ »qM*a% _ m^Nl 

Z^l^dQx) ~ L{2xf ~ Lx* ' ( ' 

i=l j—i v 1 K J 

where N m is the number of the endofullerenes with magnetic moments parallel to the magnetic 
field inside each nanotube. 

We assume that the magnetic endofullerenes inside the nanotubes form a superparamagnetic 
phase. In this case the number N m of the endofullerenes with magnetic moments parallel to the 
magnetic field is determined by the Langevin function (see, e.g., Refs. 49 and 50): 



N m = N ef C(N e fM ei ) = N, 



ef 



coth 



N ef M ef B\ k B T 



k B T J N ef M e tB 



(16) 
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where B is the magnetic field induction, 7V e f is the number of the magnetic endofullerenes inside 
the nanotube, and N e fM e f is the magnetic moment corresponding to the magnetic saturation. 

For homogeneous distribution of a load along the nanotube length the bending stiffness 
coefficient is given by [54] 

* = (") 

where the moment of inertia of the cross-section of the nanotube fully filled with the endofulle- 
renes is 

I r = ttR 4 /A. (18) 

Here R = 3nac-c/27r is the radius of the (n, n) nanotube and ac-c = 0.142 nm. 

The nanorelay turning on is possible if the magnetic force F m is larger than the elastostatic 
force F s at any deflection x of the nanotube tip. Thus the equalities of magnetic and elasto- 
static forces F s (xj) = F m (xi) and their derivatives dF s /dx\ x=X] = dF m /dx\ x=Xi determine 
the maximum value of the distance 2x a between the nanotubes (for defined N m ) for which the 
nanorelay turning on is possible (see Fig. 5). Substituting expressions (5) and (15) into these 
equalities, we obtain the following set of equations 



2 Lxf 2 Lxf 



Hence it follows that 



x, = -x , N m =--£-±— . (20) 

4 16M ef 

The expression for the average thermal oscillation amplitude in the considered case of ho- 
mogeneous distribution of a load along the nanotube length takes the form (see Eq. (10)) 

(xt-xrf = 0A243 3 -^f, (21) 

where the coefficient 3/8 is introduced into this equation because of the different bending stiff- 
ness coefficients in the cases of distributed and lumped load (cf. Eqs. (6) and (17)); x™ n = 
i? + 0.17nm. 

The number N et - of the endofullerenes in the fully filled (8,8) nanotube we estimate as 

N ei =^, (22) 

O e f 

where S e f « (0.42 + 0.34) nm is the distance between the centers of the adjacent spherical 
endofullerenes. Here 0.42 nm is the diameter of endofullerene Fe@C2o, and 0.34 nm is the 
distance between the surfaces of adjacent endofullerenes which we assume to be equal to the 
interlayer distance of graphite. 

From Eq. (20) at N m = N e f the maximal distance 2x™ ax between the nanotubes for the 
nanorelay turning on for a given nanotube lengths was obtained 



/ 16jVrfMg 
I 3 




(23) 



Substituting Eq. (20) and x = x* (according to Eq. (21)) into Eq. (16), the minimal mag- 
netic field induction B m \ n which is necessary for the nanorelay turning on can be numerically 
calculated. The dependence of the minimal magnetic field induction B m \ n on the nanotube 
length L for the nanorelay based on the (8,8) nanotubes operated at T = 5 K is shown in 
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Fig. 7. (a) The dependence of the minimal magnetic field induction B mm on the nanotube 
length L for the nanorelay operated at T = 5 K and based on the (8,8) carbon nanotubes fully 
filled with the magnetic endofullerenes Fe@C2o (left axis). Here N e f ^> 1 is the number of 
the endofullerenes in the fully filled (8,8) nanotube. The dependences of half of the minimal 
(2a;*) and maximal (2.t™ x ) distances between the nanotubes for which the nanorelay operation 
is possible on the nanotube length L (right axis); x™ n = 0.712 nm. (b) The dependence of 
the induction B necessary for the nanorelay operation with respect to the distance 2x a between 
the nanotubes for the nanotube length L = 1 /xm. Arrows indicate the values of x„, X™ 3 * and 
-Bmin shown in Fig. 7(a). The cross-section of the (8,8) nanotube (2i?j n 0.75 nm) filled with 
Fe@C2o is shown in the right part of the figure; here S e f « (0.42 + 0.34) nm and M e f = 8/xb- 



Fig. 7(a). For the (8,8) nanotubes there is no solution of Eqs. (16), (20) and (21) at room tem- 
perature (it follows from Eqs. (20) and (21) that the number N m of moments aligned with field 
necessary for the nanorelay turning on is larger than the number N e f of endofullerenes which 
the nanotube can contain (see Eq. (22)). 

When we use the (21,21) nanotubes the nanorelay can operate at room temperature (see 
Fig. 8). We assume that the endofullerenes are arranged in the hexagonal dense packing in the 
nanotube cross-section and in the form of seven chains along the nanotube axis (see the inset in 
Fig. 8(b)). It should be noted that for the nanotube of the large diameter the number N e f of the 
endofullerenes in the fully filled nanotube one can estimate as 



2 1 



7V ef = lilt 



6 SW 



L 

Sef' 



(24) 
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Fig. 8. (a) The dependence of the minimal magnetic field induction £? m j n on the nanotube length 
L for the nanorelay operated at T = 300 K and based on the (21,21) carbon nanotubes fully 
filled with the magnetic endofullerenes Fe@C2o (left axis). Here iV e f 1 is the number of 
the endofullerenes in the fully filled (21,21) nanotube. The dependences of half of the minimal 
(2x*) and maximal (2x™ ax ) distances between the nanotubes for which the nanorelay operation 
is possible on the nanotube length L (right axis); a;™" = 1.594 nm. (b) The dependence of the 
induction B necessary for the nanorelay operation with respect to the distance 2x Q between the 
nanotubes for the nanotube length L = 1 /iin. Arrows indicate the values of a;*, x™ ax and B m j n 
shown in Fig. 8(a). The cross-section of the (21,21) nanotube (2i?j n w 2.85 nm) filled with 
Fe@C2o is shown in the right part of the figure; here S e f ss (0.42 + 0.34) nm. 

where IntfiV] is the integer part of the number N, the area S e f occupied by one endofullerene 
we estimate as S e f « ^S^/A for the distance S e { w (0.42 + 0.34) nm between centers of the 
neighbor endofullerenes Fe@C2o; coefficient 7r\/3/6 is taken for the hexagonal dense packing 
of circles in the plain [63]. If the endofullerenes inside a nanotube are polymerized and form 
covalent bonds, then 6 e f will be smaller. Note that i?j n = R — 0.17 nm and R = 3nac-c/27r 
is the radius of the {n, n) nanotube. Formula (24) gives the asymptotic upper bound of N e f for 
■kRI /S* » 1. 

In Ref. 24 the existence of the magnetic endofullerene (Ho3N)@Cbo with the total magnetic 
moment of the molecule M e f = 21/xb was shown experimentally. Our calculations of opera- 
tional characteristics of the nanorelay made of (21,21) nanotubes filled with (Ho3N)@Cso at 
room temperature are presented in Fig. 9. We assume that the endofullerenes form three chains 
along the nanotube axis (see inset in Fig. 9(b)). The distance between centers of the neighbor 
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Fig. 9. (a) The dependence of the minimal magnetic field induction B m j n on the nanotube length 
L for the nanorelay operated at T = 300 K and based on the (21,21) carbon nanotubes fully 
filled with the magnetic endofullerenes (Ho3N)@Cgo (left axis). Here N e f 1 is the number of 
the endofullerenes in the fully filled (21,21) nanotube. The dependences of half of the minimal 
(2x*) and maximal (2x™ ax ) distances between the nanotubes for which the nanorelay operation 
is possible on the nanotube length L (right axis); a;™" = 1.594 nm. (b) The dependence of the 
induction B necessary for the nanorelay operation with respect to the distance 2x Q between the 
nanotubes for the nanotube length L = 0.75 /im. Arrows indicate the values of x„, x™ ax and 
-Bmin shown in Fig. 9(a). The cross-section of the (21,21) nanotube (2i?j n « 2.85 nm) filled 
with (Hc>3N)@C8o is shown in the right part of the figure; here <5 e f ~ (0.8 + 0.34) nm and 
M ef = 21^ B - 

endofullerenes (Ho3N)@Cgo is <5 e f ~ (0.8 + 0.34) nm. Here 0.8 nm is the diameter of endofulle- 
rene (Ho3N)@Qo, and 0.34 nm is the distance between the surfaces of adjacent endofullerenes 
which we assume to be equal to the interlayer distance of graphite. 

5 CONCLUSION 

In this paper, we applied the DFT to calculate the structural and energy characteristics of the 
smallest magnetic endofullerene Fe@C2o- The ground state of Fe@C2o was found to be the 
septet state, and the magnetic moment of Fe@C2o was estimated to be 8/ib- The characteristics 
of the equilibrium position and motion of this endofullerene inside the (8,8) nanotube were 
studied in the framework of the semiempirical approach. The calculated barriers for the motion 
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and rotation show that the endofullerene freely moves and rotates inside the nanotube at room 
temperature. 

The scheme of the magnetic nanorelay based on cantilevered nanotubes filled with magnetic 
endofullerenes is elaborated. The proposed nanorelay is turned on as a result of bending of the 
nanotubes by a magnetic force. Operational characteristics of the nanorelay based on the (8,8) 
and (21,21) nanotubes fully filled with the proposed here endofullerenes Fe@C2o and on (21,21) 
nanotube fully filled with observed (Ho3N)@Cgo endofullerenes are calculated. It is shown that 
the nanorelay based on the (21,21) nanotubes fully filled with the endofullerenes Fe@C2o or 
(Ho3N)@Cgo can operate at room temperature for the nanotubes of about 1 fim in length. 

The proposed nanorelay can be used also for the measurements of magnetic moments of 
endofullerenes or other magnetic molecules inside single-walled carbon nanotubes and for in- 
vestigation of elastic properties of nanotubes. 
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